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ABSTRACT 

This is the third of a series of papers in which we derive simultaneous constraints on 
cosmological parameters and X-ray scaling relations using observations of the growth 
of massive, X-ray flux-selected galaxy clusters. Our data set consists of 238 clusters 
drawn from the R OS A T All-Sky Survey, and incorporates extensive follow-up obser- 
vations using the Chandra X-ray Observatory. Here we present improved constraints 
on departures from General Relativity (GR) on cosmological scales, using the growth 
index, 7, to parameterize the linear growth rate of cosmic structure. Using the method 
of Mantz et al. (2009a), we simultaneously and sclf-consistently model the growth of 
X-ray luminous clusters and their observable-mass scaling relations, accounting for sur- 
vey biases, parameter degeneracies and the impact of systematic uncertainties. Such 
analysis of the survey and follow-up data is crucial, else spurious constraints may be 
obtained. We combine the X-ray cluster growth data with cluster gas mass fraction, 
type la supernova, baryon acoustic oscillation and cosmic microwave background data. 
We find that the combination of these data leads to a tight correlation between 7 and 
the normalization of the matter power spectrum, og ■ Consistency with GR requires a 
measured growth index of 7 ~ 0.55. Under the assumption of self-similar evolution and 
constant scatter in the cluster observable-mass scaling relations, and for a spatially 
flat model with a cosmological constant, we measure 7(crs/0.8) 6 8 = 0.55^0 ioi with 
allowed values for erg in the range 0.79 to 0.89 (68.3 per cent confidence limits). Relax- 
ing the assumptions on the scaling relations by introducing two additional parameters 
to model possible evolution in the normalization and scatter of the luminosity-mass 
relation, we obtain consistent constraints on 7 that are only ~ 20 per cent weaker 
than those above. Allowing the dark energy equation of state, w, to take any constant 
value, we simultaneously constrain the growth and expansion histories, and find no 
evidence for departures from either GR or the cosmological constant plus cold dark 
matter paradigm. Our results represent the most robust consistency test of General 
Relativity on cosmological scales to date. 
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1 INTRODUCTION 

Since Riess et al. (1998) and Perlmutter et al. (1999) first 
demonstrated, using observations of type la supernovae 
(SNIa), that the expansion history of the Universe is ac- 
celerating at late times, improved SNIa studies have con- 
tinuously corroborated these results (Garnavich et al. 1998; 
Schmidt et al. 1998; Knop et al. 2003; Tonry et al. 2003; 
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Barris et al. 2004; Riess et al. 2004, 2007; Astier et al. 2006; 
Wood-Vasey et al. 2007; Kowalski et al. 2008; Hicken et al. 
2009). At the same time, measurements of anisotropies in 
the cosmic microwave background (CMB) with the Wilkin- 
son Microwave Anisotropy Probe (WMAP Spergel et al. 
2003, 2007; Komatsu et al. 2009; Dunkley et al. 2009, 
and companion papers) and other CMB experiments 
(Readhead et al. 2004; Jones et al. 2006; Reichardt et al. 
2009; Chiang et al. 2010; Gupta et al. 2009) have tightly 
constrained other key cosmological parameters. A variety 
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of other cosmological data sets have also been used to in- 
dependently measure cosmic acceleration. Measurements of 
the gas mass fraction (/gas) in galaxy clusters, having ear- 
lier shown that the mean matter density of the Universe is 
low, Q m ~ 0.25 (see e.g. White et al. 1993), have directly 
confirmed the effects of cosmic acceleration, at comparable 
significance to that seen in SNIa data (see e.g. Allen et al. 
2004, 2008). Measurements of baryon acoustic oscillations 
(BAO) in galaxy surveys (see e.g. Eisenstein et al. 2005; 
Percival et al. 2007, 2010), the correlation of the low mul- 
tipoles of the CMB with large scale structure observa- 
tions from various surveys (see e.g. Scranton et al. 2003; 
Fosalba et al. 2003; Cabre et al. 2006; Giannantonio et al. 
2008; Ho et al. 2008), and weak gravitational lensing by 
large scale structure (Schrabback et al. 2009) have also 
shown evidence for cosmic acceleration. The combination 
of subsets and/or all of these data have led to the establish- 
ment of the ACDM paradigm, in which the energy density of 
the Universe is currently composed of ~ 5 per cent baryons, 
~ 20 per cent cold dark matter (CDM) and ~ 75 per cent 
dark energy, the last of which drives cosmic acceleration and 
has identical characteristics to Einstein's cosmological con- 
stant, A. 

Together, the above experiments robustly show that the 
Universe is accelerating. However, the ability of such data 
to probe the underlying cause of this acceleration is limited. 
All of the above constraints on dark energy are primarily 
driven by its effects on the background geometry. Such data 
alone cannot distinguish acceleration due to a true dark en- 
ergy component with negative pressure from modifications 
to the standard theory of gravity, i.e. Einstein's General Rel- 
ativity (GR). Recently, however, Mantz et al. (2008, 2009a, 
hereafter Paper I), and Vikhlinin et al. (2009) have pre- 
sented new constraints on dark energy from measurements 
of the growth of cosmic structure, as evidenced in X-ray 
flux-selected cluster samples. These experiments are sensi- 
tive not only to the evolution of the mean energy density 
background, but also to the evolution of the density per- 
turbations with respect to this background. This work has 
provided clear, independent confirmation of the effects of 
dark energy in slowing the growth of X-ray luminous galaxy 
clusters. 

The observed evolution of density perturbations pro- 
vides a sensitive probe of the underlying theory of gravity 
and the clustering 1 properties of dark energy. Rapetti et al. 
(2009, hereafter R09) exploited this sensitivity to constrain 
departures from the predicted cosmic growth rate for GR 
using the convenient parameterization f2 m (z) 7 , for which 
GR has 7 ~ 0.55 (Peebles 1980; Wang & Steinhardt 1998; 
Linder 2005; Huterer & Linder 2007; Di Porto & Amendola 
2008; Sapone & Amendola 2007; Gannouji & Polarski 2008; 
Wei 2008; Linder & Cahn 2007; Nesseris & Perivolaropoulos 
2008). Using the growth of structure analysis of Mantz et al. 
(2008), R09 found no evidence for deviations from either 
GR or ACDM. Recent results from the combination of 



1 Note that for models other than the cosmological constant, dark 
energy is expected to couple with gravity, i.e. cluster (see e.g. 
Hu 2005; Mota et al. 2007). For specific dark energy and modi- 
fied gravity models, see the reviews of Copeland et al. (2006) and 
Frieman et al. (2008). 



other cosmological data sets are also consistent with GR 
(Daniel et al. 2010; Zhao et al. 2010; Reyes et al. 2010). 

This is the third of a series of papers in which, for the 
first time, we employ a fully self-consistent analysis of the 
growth of massive clusters that combines current X-ray clus- 
ter surveys and deep, pointed, follow-up observations to si- 
multaneously constrain both cosmological parameters and 
observable-mass (luminosity-mass and temperature-mass) 
scaling relations. Importantly, the improved method, which 
is described in Paper I, properly accounts for all selection 
biases, covariances and parameter degeneracies, and mod- 
els fully the impact of systematic uncertainties. Here, we 
utilize this improved method (hereafter XLF, as defined in 
Paper I) to re-investigate the simultaneous constraints that 
can be placed on 7 and the background evolution (expansion 
history). For the background, we use two reference expan- 
sion models: flat ACDM, parameterized by the mean matter 
density, fl m ; and flat iuCDM, parameterized by $l m and a 
constant dark energy equation of state, w. Our constraints 
on 7 arise primarily from the XLF experiment, which uses 
the cluster samples of Ebeling et al. (1998, ROSAT Bright- 
est Cluster Sample; BCS), Bohringer et al. (2004, ROSAT- 
ESO Flux Limited X-ray sample; REFLEX), Ebeling et al. 
(2001, 2009 in preparation, Bright MAssive Cluster Survey; 
Bright MACS), plus extensive X-ray follow-up data from 
the Chandra X-ray Observatory and ROSAT (for details see 
Mantz et al. 2009b, hereafter Paper II), although we also 
utilize the (currently) small additional constraining power 
available from measurements of the Integrated Sachs- Wolfe 
(ISW) effect at low multipoles in the CMB. 

We emphasize that the models used for the evolution 
of the background and density perturbations are purely 
phenomenological, encompassing both modified gravity and 
clustering dark energy as possible sources for cosmic ac- 
celeration. Our analysis allows for a simple, but elegant 
test for departures from the standard GR and ACDM 
paradigms. Using the combination of XLF, / ga s, SNIa, BAO 
and CMB, we demonstrate good agreement with the stan- 
dard GR+ACDM model and provide tight constraints on 
the model parameters. We show that our results are robust 
against reasonable assumptions regarding the evolution of 
galaxy cluster scaling relations. We stress that in order to 
obtain robust results from current and future XLF studies, 
it is essential to employ a consistent analysis method, such 
as the one used here (see details in Paper I), in order to 
properly account for selection biases, covariances, param- 
eter degeneracies and systematic uncertainties. Otherwise, 
spuriously tight constraints may be obtained. 



2 COSMOLOGICAL MODEL 
2.1 Linear growth rate 

Galaxy clusters are the largest and rarest virialized objects 
in the Universe. Their abundance as a function of redshift 
provides an extremely sensitive probe of the underlying cos- 
mology. To predict the cumulative number density of galaxy 
clusters, n(M,z), at a given mass, M, and redshift, z. we 
calculate the evolution of density fluctuations using linear 
perturbation theory, accounting for non-linear effects using 
results from N-body simulations. 
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The cumulative mass function of dark matter halos can 
be calculated as 



n(M,z) = 



, . p m diner 1 , 



(1) 



where p m is the mean comoving matter density, a is the 
variance of the linearly evolved density field (as defined 
in Equation 5), and f(a) is a fitting formula obtained 
from theory or N-body simulations (Press & Schechter 1974; 
Bond et al. 1991; Sheth & Tormen 1999; Jenkins et al. 
2001; Evrard et al. 2002). Remarkably, f{o) encompasses 
linear and non-linear effects in such a way that its form is 
approximately independent of the cosmology assumed (see 
Section 2.3). 

Modifications to GR at large scales have been proposed 
as alternatives to dark energy to explain cosmic accelera- 
tion at late times. In general, in addition to changing the 
expansion history, such modifications also affect the growth 
of cosmic structure. As in R09, we parameterize the growth 
rate, f(a), using (Peebles 1980; Wang & Steinhardt 1998; 
Linder & Cahn 2007) 
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where f2 m (a) = S7 m a~ 3 /Z?(a) 2 . Here, Q m is the present, 
mean matter density in units of the critical density, E(a) = 
H(a)/Ho the evolution parameter, H(a) the Hubble param- 
eter (with Ho its present-day value), and 7 the growth index. 

We calculate the evolution of the linear matter density 
contrast, 5 = 5p m /p m , i.e. the ratio of matter density fluc- 
tuations to the cosmic mean, by solving Equation 2. As an 
initial condition for this differential equation, we match the 
value of S at early times (z t = 30) to that of GR. Given 5(z), 
we obtain the growth factor between z% and z, 



D(z) = 



S(z t ) 



(3) 



Note that z t = 30 is well within the dark matter dominated 
era, when the expansion is expected to be decelerating. At 
this time f(a) is close to 1, independent of 7. Using D(z), 
we calculate the linear matter power spectrum for a given 
wavenumber, k, and redshift, 



P(k,z) oc k Us T 2 {k,z t )D{zf 



(4) 



where T(k,z t ) is the matter transfer function of GR at red- 
shift z t , and n s the primordial scalar spectral index. Given 
P(k,z), the variance of the linearly evolved density field, 
smoothed by a spherical top-hat window function of comov- 
ing radius R enclosing mass M = 47rp m J? 3 /3, is then 

a 2 (M,z) = — 5- / k 2 P(k,z)\W M (k)\ 2 dk, (5) 
StH Jo 

where Wyi(k) is the Fourier transform of the window func- 
tion. Using this value for a(M,z) in Equation 1, we obtain 
n(M, z). 



2.2 Background evolution 

Equation 2 depends on the evolution of the background en- 
ergy density, through the parameter E(a). Assuming an ex- 
pansion history appropriate for ACDM, it has been shown 
(see e.g. Linder & Cahn 2007; Nesseris & Perivolaropoulos 



2008) that, for 7 ~ 0.55, the growth rate calculated using 
Equation 2 corresponds remarkably accurately to that pre- 
dicted by GR. Linder & Cahn (2007) also showed that 7 and 
w are only weakly correlated (see also R09) and, therefore, 
that one can measure the cosmic growth and expansion rates 
almost independently. These convenient properties of the 7- 
model provide us with a simple but powerful consistency 
test for GR. 

Here, we perform this test using the new cluster XLF 
data (Papers I and II) in combination with / gas , SNIa, BAO 
and CMB constraints. The latter four data sets primarily 
constrain the background evolution of the Universe, i.e. the 
expansion history. As a background expansion model, we 

„„„2 



E(a)= [fi m oT 3 + (l 
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where w is the dark energy equation of state parameter and 
w = — 1 corresponds to flat ACDM. Note that we use this 
model only as a phenomenological description of the expan- 
sion history. We do not assume the existence of dark energy. 
Our model allows us to simultaneously test for departures 
from ACDM in the expansion history (w / —1) and depar- 
tures from GR in the growth history (7 7^ 0.55). Deviations 
from these benchmarks can be caused either by clustering 
dark energy or modified gravity. 



2.3 Non-linear model 

The fitting formula f(a) from Jenkins et al. (2001) has 
been tested for a wide range of background dark energy 
cosmologies and has been shown to provide a "universal" 
description, within ~ 20 per cent uncertainty (see e.g. 
Kuhlen et al. 2005; Lokas et al. 2004; Klypin et al. 2003; 
Linder & Jenkins 2003; Mainini et al. 2003; Maccio et al. 
2004; Francis et al. 2009; Grossi & Springel 2009) 3 . How- 
ever, in these tests, typically only the background evolu- 
tion was modified, while the linear density perturbations 
of the new component (dark energy) were neglected. Re- 
cently, Jennings et al. (2010) have presented results that ac- 
count for the effects of the density fluctuations provided by 
quintessence models on the linear matter power spectrum, 
finding only a small effect on f(a) for the masses and red- 
shifts relevant here. 



2 R09 examined three expansion models: flat ACDM, flat miCDM 
and non-flat ACDM. For the latter, they found that there is neg- 
ligible covariance between fl^, the curvature energy density, and 
7. Since this test is computationally expensive, we do not repeat 
it here. 

3 Schmidt et al. (2009) explored an f(R) gravity model and 
showed that for values of their model parameter compatible with 
Solar Systems tests (/ro ~ 10~ 6 ), the halo mass function is 
in agreement with that of GR for very massive clusters (within 
the quoted level of uncertainty of the simulations). However, for 
more extreme values of the model parameter and for less mas- 
sive clusters (for which chameleon effects arc important), sig- 
nificant differences can arise. For the self-accelerating branch of 
the braneworld gravity model DGP (Dvali et al. 2000), Schmidt 
(2009a) and Chan & Scoccimarro (2009) predict a significant sup- 
pression of the halo mass function relative to GR, especially for 
massive objects. Schmidt (2009b) predicts an enhancement for 
the normal branch of this model. 
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The present study aims to test the consistency of cur- 
rent observations with GR (and ACDM). Therefore, we use 
an f(a) corresponding to GR, seeking to determine whether 
values of 7 ~ 0.55 (and w = — 1) are preferred by the data. 
Beyond the mass function of Jenkins et al. (2001), several 
authors have investigated the inclusion of an explicit depen- 
dence on redshift in the fitting formula f(a, z) (Lukic et al. 
2007; Reed et al. 2007; Cohn & White 2008). Following Pa- 
per I, we use the f(a,z) function of Tinker et al. (2008), 
determined from a large suite of simulations, 



A 



+ 1 



(-) 



Here, parameters have the redshift dependence x(z) = 
a:o(l + z) £0!x , with x £ {A, a, b, c}. The values for xo and cv x 
are given by Tinker et al. (2008). As described in Paper I, we 
use the parameter s to account for residual systematic uncer- 
tainties in the evolution of f(a, z). We marginalize over the 
uncertainties in the parameters of Equation 7, accounting 
for their covariance. We also account for additional system- 
atic uncertainties due to e.g. the presence of baryons (for 
details, see Paper I). Note, however, that, as we show in Pa- 
per I, the statistical and systematic uncertainties in f(a, z) 
are subdominant in the analysis. We have verified that e is 
essentially uncorrelated with 7. 



EVOLUTION OF THE SCALING 
RELATIONS 



In this paper we present results both for pure self- 
similarity and constant scatter (/3| m = 0, a' lm = 0), and 
also allowing for additional evolution in the luminosity-mass 
relation and its scatter, through f3 2 m and a' tm . As we shall 
show, our improved analysis method (of Paper I) and the 
inclusion of high quality follow-up data for a significant frac- 
tion of clusters allows us to address even general questions of 
this type. We emphasize that for such analysis it is essential 
to simultaneously model the mass function, scaling relations, 
growth history and the impact of systematic uncertainties 
fully, else spurious constraints may be obtained. 

The data and analysis of Paper I and II also in- 
clude measured temperatures to simultaneously constrain 
the temperature-mass relation. In Paper II, we show that a 
simple power law 



(t(m))=^ m + p[ m m, 



(10) 



where t = log 10 (fcTsoo/keV), without additional evolution 
parameters, such as /?2 m and a' tm (which are defined equiv- 
alently to the corresponding parameters in the luminosity- 
mass relation, /3 2 m and a' tm ), is sufficient to describe the 
data. Note also that, for clusters with kT > 3keV, such as 
those used here, the flux within the 0.1 — 2.4 keV band is 
nearly independent of the temperature, and therefore fi^ 1 
and a' tm are essentially uncorrelated with 7. 5 Thus, we do 
not further consider these parameters, and assume that the 
temperature-mass relation evolves self-similarly, with a con- 
stant log-normal scatter 6 . 



To compare the XLF data with the theoretical predictions of 
the mass function, we need to relate mass, M, to our observ- 
ables: the X-ray luminosity, L, and average temperature, T, 
of the clusters. In Paper II, we describe the follow-up X-ray 
observations used to determine these relations. As discussed 
in Paper I, we simultaneously and self-consistently constrain 
both cosmological (see Section 2) and scaling relation pa- 
rameters. 

Following the notation of Paper II, we model the evo- 
lution of the luminosity-mass scaling relation as 



(l(m)) = ^ + p[ m m + Pi m log 10 (l + z) , 



(8) 



with a log-normal, and possibly evolving, intrinsic scatter of 
the luminosity at a given mass 4 



crtm{z) = o lm {\ + a' tm z) 



(9) 



Here, I = log 10 [L 50 oi?(z)~ 1 /10 44 ergs -1 ] and m = 
log 10 [-E(z)M 50 o/10 15 M e ], where the subscript 500 refers to 
quantities measured within radius rsoo, at which the mean, 
enclosed density is 500 times the critical density of the Uni- 
verse at redshift z. The factors of E(z) are required to ac- 
count for the background evolution of the critical density. 
Fixing /3| m = 0, one has "self-similar" evolution of the scal- 
ing relation (Kaiser 1986; Bryan & Norman 1998), deter- 
mined entirely by the E{z) factors. Fixing a' lm = 0, one 
has constant scatter. In Paper II, we show that departures 
from self-similar evolution and evolution in tje m (z) are not 
required by current data. 



4 In the notation of R09, /3^ m , /3f m , /3| m , <r tn 
spond to log 1 Q(Afo), P, Ci Vo and rj z , respectively. 



and <r' e corre- 



4 ISW EFFECT 

Our analysis also accounts for the small, but non-negligible 
additional constraining power on 7 that is available from 
measurements from the ISW effect. The low multipoles of 
the CMB are sensitive to the growth of structure through 
this effect. The time- varying gravitational potentials of large 
scale structures contribute a net effect on the energy of the 
photons crossing them. For these photons, we calculate their 
contribution to the temperature anisotropy power spectrum 
as (Weller & Lewis 2003) 



Af W (fc) = 2 



dt e 



4> ji [k{t - t )] 



(11) 



where t is the conformal time, r the optical depth to reion- 
ization, ji (x) the spherical Bessel function for the multipole 
I, and <j)' the conformal time variation of the gravitational 
potential. The latter can be calculated in terms of 7 using 
the Poisson equation, k 2 (f> = —AirGa 2 5p m , as (R09) 



4nGa 2 
k 2 



H5p m [1 - O m (a) 



(12) 



where H is the conformal Hubble parameter. Since we 
are performing a consistency test of GR, we assume that, 
as in GR, the contributions of the anisotropic stress and 
energy flux are negligible (for details on these terms see 



5 We have explicitly verified this to be the case. 

6 Our analysis includes a correlation parameter between the in- 
trinsic scatters of the luminosity-mass and temperature-mass re- 
lations, Pit m ( see details in Paper II). 
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Figure 1. 68.3 and 95.4 per cent confidence contours in the n m ,7 (left panel) and erg, 7 (right panel) planes for an assumed flat 
ACDM background model, using the combination of XLF, /gas, SNIa, BAO and CMB data. The gold, smaller contours assume self- 
similar evolution of the observable-mass scaling relations and constant scatter (/3| m = 0; <?' trn = 0). The blue, larger contours allow for 
departures from self-similarity and redshift evolution in the scatter of the luminosity- mass relation (Section 3; Paper II). The horizontal, 
dashed lines mark 7 = 0.55, the growth index for GR. 



Challinor & Lasenby 1999) 7 . As an initial condition, at z = 
2 we match the Aj sw (fc) to that of GR, using Equation 12 
to evolve the model to z = 0. Unless stated, for results in- 
cluding CMB data, we include the constraints on 7 available 
from the ISW effect. 



5 DATA ANALYSIS 

To constrain the cosmic growth rate we use three wide- 
area cluster samples drawn from the ROSAT all-sky sur- 
vey: BCS (z < 0.3; northern sky; F x (0.1 - 2.4 keV) > 4.4 x 
10" 12 ergs" 1 cm" 2 ), REFLEX (z < 0.3; southern sky; F x > 
3.0 x 10" 12 ergs" 1 cm" 2 ), and Bright MACS (0.3 < z < 0.5; 
~ 55 per cent sky coverage; fx > 2 x 10~ 12 ergs~ 1 cm" 2 ). 
In order to keep systematic uncertainties to a minimum, 
for all three surveys, we impose a lower luminosity cut of 
2.5 x 10 44 /if 2 ergs" 1 (0.1 - 2.4 keV). For the BCS, this gives 
a total of 78 clusters; for REFLEX, 126 clusters; and for 
Bright MACS, 34 clusters. In total we have 238 clusters. 

To constrain the background evolution we also use mea- 
surements of the gas mass fraction (/gas) in 42 massive, 
dynamically relaxed clusters from Allen et al. (2008), as 
well as the compilation of 307 SNIa from Kowalski et al. 
(2008), and BAO data (Percival et al. 2007; Colless et al. 
2001, 2003; Adelman-McCarthy et al. 2007). The latter con- 
strain the ratio of the sound horizon to the distance scale 
at z = 0.25 and 2 = 0.35. We also use the five-year WMAP 
data (Dunkley et al. 2009; Komatsu et al. 2009, and com- 
panion papers). 

7 For particular alternative gravity models such as DGP, these 
terms will not be negligible (Hu & Sawicki 2007) and the Poisson 
equation will need to be modified (Hu 2008). 



We calculate the posterior probability distributions of 
all parameters using the Metropolis Markov Chain Monte 
Carlo (MCMC) algorithm, as implemented in the Cos- 
MOMC 8 code of Lewis & Bridle (2002). We use a modified 
version of this code that includes additional modules to cal- 
culate the likelihood for our two cluster experiments: / ga s 9 
(Rapetti et al. 2005; Allen et al. 2008) and the XLF (Pa- 
per I). For analyses without CMB data, we use standard 
Gaussian priors on both the present-day Hubble constant, 
H = 72 ± 8 km s^Mpc" 1 (Freedman et al. 2001), and 
the mean physical baryon density, Q h h 2 = 0.0214 ± 0.0020 
(Kirkman et al. 2003), from Big Bang Nucleosynthesis stud- 
ies. When including the CMB data, we do not use these 
priors. 

When using CMB data, instead of Ho we fit 6. The 
latter is the (approximate) ratio of the sound horizon at last 
scattering to the angular diameter distance, which, as shown 
by Kosowsky et al. (2002), is less correlated than Hq with 
other parameters. In all analyses, in addition to Hq or 6, 
and Qhh , we fit for the mean physical dark matter density, 
Q c h 2 ; the logarithm of the adiabatic scalar amplitude, ln(A s ) 
(which is related to ag) 10 ; the growth index, 7; and, when 
stated, the kinematical w parameter 11 . For the analyses with 
CMB data, we also fit the optical depth to reionization, r, 
the adiabatic scalar spectral index, n s , and marginalize over 
the amplitude of the Sunyaev-Zel'dovich effect from galaxy 
clusters, < Asz < 2 (Spergel et al. 2007). For analyses 

8 http:/ /cosmologist.info/cosmomc/ 

9 http:/ /www. stanford.edu/~drapetti/fgas_module/ 

10 <t| is the 2 = variance in the density field at scales of 
8/1" 1 Mpc (see equation 5). 

11 Recall that we use w purely as a kinematic parameter, allowing 
us to model the cosmic expansion history conveniently. 
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Figure 2. 68.3 and 95.4 per cent confidence contours in the mi, 7 plane for the flat uiCDM background expansion model, using the 
combination of XLF, / ga s, SNIa, BAO and CMB data. The results assume self-similar evolution of the scaling relations and constant 
scatter (/3| m = 0; = 0). The horizontal, dashed line marks 7 = 0.55, the growth history for GR. The vertical, dotted-dashcd line 
marks w = —1, the expansion history for ACDM. The results are simultaneously consistent with GR and ACDM. 



without CMB data, we fix n s to 0.95 since, for such analyses, 
ris is degenerate with as (see Paper I). 

6 SIMULTANEOUS CONSTRAINTS ON THE 
GROWTH AND EXPANSION HISTORIES 

6.1 Results for the ACDM expansion model 

Figure 1 shows the joint constraints in the fi m , 7 (left panel) 
and (78, 7 (right panel) planes for a fiat ACDM expansion his- 
tory. The gold (smaller) contours show the constraints for 
the minimal, self-similar and constant-scatter scaling rela- 
tion model (/3 l 2 m = and a' lm = 0). 

Figure 1 shows that, for the ACDM background model, 
Qm is not strongly correlated with 7. This supports the idea 
that 7 is a useful parameterization to measure departures 
from the growth rate of GR, approximately independent of 
assumptions regarding the background evolution, parame- 
terized in this case by f2 m . 

The right panel of Figure 1 shows a clear correlation 
between as and 7, with a correlation coefficient of p = 
-0.87. From the MCMC samples we measure 7(08 /0.8) 6,8 = 
0.55lo jq, with allowed values for as in the range 0.79 < 



as < 0.89 (marginalized 68.3 per cent confidence limits) 12 . 
For as = 0.8, we measure 7 = 0.55 (GR; marked by a hor- 
izontal, dashed line in Figure 1). Our results are consistent 
with GR at the 68.3 per cent confidence level. The observed 
correlation between erg and 7 shows that the addition of 
precise, independent knowledge of as should significantly 
improve the precision of the constraints on 7 (see also Sec- 
tion 6.4). 

The marginalized constraints on a single interesting pa- 
rameter are summarized in Table 1. 



6.2 Results for the uiCDM expansion model 

Figure 2 shows the joint constraints in the w, 7 plane for the 
toCDM background cosmology (w constant), assuming self- 
similar evolution and constant scatter (/?2 ?n = 0; a' lm = 0). 
The horizontal, dashed and vertical, dotted-dashed lines in- 
dicate 7 = 0.55 (GR) and w = -1 (ACDM), respectively. 
Importantly, our results are simultaneously consistent with 



We calculate the exponent 6.8 finding the best constrained 
eigenvector for the estimated covariance matrix between In 7 and 
In erg . 



© 2008 RAS, MNRAS 000, 1-11 



Testing General Relativity with the observed growth of massive clusters 7 




Figure 3. 68.3 and 95.4 per cent confidence contours in the /?2 >7 (l e ft panel) and CT< m ,7 (right panel) planes for the flat ACDM 
background evolution model, from the combination of XLF, /gas, SNIa, BAO and CMB data. The green, smaller contours show the 
constraints obtained with either /3j m (left panel) or u' lm (right panel) allowed to be a free parameter. The blue, larger contours show 
the results when both parameters, /3| m an d a 't m i are allowed to be free. The horizontal, dashed lines mark 7 = 0.55 (GR). The vertical, 
dottcd-dashed lines mark the self-similar (/3| m = 0; left panel) and constant scatter (c' em = 0; right panel) conditions. 



both GR and ACDM at the 68.3 per cent confidence level 
(see Figure 2). For this expansion model, we again observe 
a correlation between as and 7, although somewhat less 
(p = —0.69) than for the ACDM case. The marginalized 
constraints on a single interesting parameters are summa- 
rized in Table 1. 



6.3 Testing for departures from self-similarity 

The blue (larger) contours in Figure 1 show the constraints 
obtained for the ACDM background model when allowing 
/?2 m and a' tm to be free. Remarkably, using this (signifi- 
cantly) more complex model for the luminosity-mass rela- 
tion, we obtain constraints on 7 that are only ~ 20 per cent 
weaker than those for the minimal, self-similar, constant- 
scatter model (gold, smaller contours). The primary reason 
for this robustness is the comprehensive nature of the follow- 
up data in the XLF experiment, which tightly constrain the 
scaling relations and their evolution (Paper II). 

Figure 3 shows the weak correlations between 7 and 
the additional evolution parameters, /3| m (left panel) and 
a 'im (right panel). The blue (larger) contours in both panels 
are the constraints obtained with both parameters free; the 
green (smaller) contours show the results for models in which 
only a single additional parameter (either /3| m or a' im ) is 
allowed to vary (fixing the other to 0; see also Table 1). 

In Paper II, we show that, for GR models and a 
flat ACDM background, current X-ray data at 2 < 0.5 
(survey+follow-up) do not require additional evolution of 
the scaling relations beyond the self-similar, constant scat- 
ter model (as implemented using the parameters /3 2 m , 02™, 
a' lm and a' tm ; see Section 3). Using the Deviance Informa- 
tion Criterion (DIC) of Spiegelhalter et al. (2002) (see defi- 




-1.0 -0.5 0.0 0.5 1.0 1.5 
7 

Figure 4. Marginalized constraints on 7 for the flat ACDM back- 
ground expansion model, from the XLF data alone (green, dottcd- 
dashed line), from the combination of SNIa, / ga s, BAO and CMB 
data (red, dashed line), and from all the data sets combined, in- 
cluding the XLF (blue, solid line). It is clear that the XLF domi- 
nates the constraints on 7 and that, alone, it provides significantly 
tighter constraints than those from the combination of the other 
data sets. 



nition in Paper II) we have verified that the minimal model 
{fi\ m = 0, <j'i m — 0) remains an adequate description of the 
data, when 7 is included as a parameter in the analysis. 
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Figure 5. 68.3 and 95.4 per cent confidence contours in the f2 m ,7 (left panel) and eg, 7 (right panel) planes for the flat ACDM 
background expansion with self-similar evolution and constant scatter in the scaling relations. Results are shown for the following 
combinations of data sets: XLF+/ gas (red contours), XLF+/ gas +SNIa (green contours), XLF+/ gas +SNIa+BAO (blue contours), and 
XLF+/gas+SNIa+BAO+CMB (gold, smallest contours). This figure does not include the additional constraints on 7 from the ISW 
effect. The horizontal, dashed lines mark 7 = 0.55 (GR). 



Table 1. Marginalized 68.3 per cent confidence constraints for various subsets of the data, for the flat ACDM (w = —1) and flat uiCDM 
(w constant) background expansion models. Note: "these results do not include the constraints on 7 from the ISW effect. 



Data 




^ 2 m 










w 


7 


XLF+/gas 





91 ,+0.036 
z± -0.041 





oc+0.07 
oo -0.06 








-1 


42+ 20 
u '^ z -0.16 


XLF+/ g a S +SNIa 





260+ 040 
ZDU -0.025 





oq+0.06 
o,3 -0.07 








-1 


n 4 n+0.25 
U ' 4J -0.15 


XLF+/ ga8 +SNIa+BAO 





9zl7 +0.028 
-0.024 





oc+0.05 
oo -0.06 








-1 




XLF+/ gas +SNIa+CMB 





ZDD -0.021 





ar;+0.04 
°°-0.06 








-1 


o-4o±g:« 


XLF+/ gas +SNIa+BAO+CMB a 





24 o+0.017 
z ^ -0.012 





84+0-06 

°^-0.06 








-1 


42+0- 15 

u '^ z -0.15 


XLF+/ gas +SNIa+BAO+CMB 





249 +0.015 
zy,a -0.013 





84+0-05 

°^-0.05 








-1 


u '^ u -0.12 


55 





9 r- 9 +0.015 
zoz -0.014 





o 9 +0.05 
oz -0.05 


n 16+ - 44 
U -' ;>D -0.51 





-1 


n 44+0- 17 
u -^-0.13 







949 +0.018 
z * -0.010 





o 9 +0.06 
-0.04 





-0.52±g;g 


-1 


45+ 14 
u '^°-0.13 


55 





24 q+0.017 
Z4y -0.011 





84 -0.06 


n S3+ 72 


-o.82t°; 5 37 


-1 


n 4n+ 017 
U - 4U -0.13 


55 





nn +0.016 
zo± -0.014 





QQ + 0.06 

° a -0.04, 








-0.98«;« 7 


39+°- 14 
u - oa -0.13 



6.4 The impacts of the different data sets 

The green, dotted-dashed line in Figure 4 shows the 
marginalized constraints on 7 from the XLF data alone 13 . 



13 Note that the XLF data include / gas measurements for the six 
lowest redshift clusters in the sample of Allen et al. (2008), which 
serve to normalize the M gas mass proxy (see details in Paper I). 
For the XLF data alone, we obtain the same constraints on 7 and 
(78 than those from the XLF+/ gas (full sample) data (see Tabic 
1), and only ~ 10 weaker constraints on f2 m . 



The red, dashed line shows the constraints from the com- 
bination of SNIa, / gas , BAO, and CMB data, excluding the 
XLF. It is readily apparent that the constraints on 7 from 
the combination of all other data sets, excluding the XLF, 
are significantly weaker than those from the XLF data alone. 
Without the XLF data, the constraints on 7 are primarly 
provided by the ISW effect, and we obtain 7 = 0.38lo'g|. 
Combining the XLF and other data, we obtain the blue, solid 
curve and constraints on 7 more than three times tighter 
than without the XLF (see Table 1). 

Although the XLF data dominate the constraints on 7, 
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the other data (/gas, SNIa, BAO, and CMB) provide com- 
plementary information on other cosmological parameters. 
These additional constraints break degeneracies with 7 and 
help to exploit fully the ability of the XLF data to con- 
strain 7. To demonstrate the impacts of these data sets, 
Figure 5 shows the constraints in the Q m ,7 (left panel) and 
erg, 7 (right panel) planes for various subsets of data. The 
red contours show the constraints from the XLF+/ gas data 
(i.e. only cluster data); the green contours show the results 
from adding SNIa data (i.e. XLF+/ gas +SNIa); the blue con- 
tours from adding BAO to the rest; and the gold (smallest) 
contours from adding CMB data 14 . 

The left panel of Figure 5 demonstrates again the ab- 
sence of any strong correlation between f2 m and 7 (see also 
Section 6.1). These results (see also Table 1) show that sim- 
ply improving our knowledge of fi m , and therefore the ex- 
pansion history (by e.g. adding BAO to XLF+/ gas -|-SNIa 
data), will not necessarily improve the constraints on 7. 

The right panel of Figure 5 shows the constraints in 
the crg,7 plane. We see that the combination of the CMB 
data with the other data sets considerably strengthens the 
correlation between these parameters. The reason for this is 
that the inclusion of the CMB data significantly reduces the 
range of allowed values not only for the expansion history 
but also for other relevant cosmological parameters, placing 
tight constraints on, e.g., Q.bh 2 , fl c h 2 , n a , r, Ho, and A B . 
These constraints lead to a tight constraint on a$(z) at early 
times (z ~ 1100), where GR is assumed (i.e., z > z t ). 



7 CONCLUSIONS 

We have used measurements of the evolution of massive 
galaxy clusters from wide-area X-ray cluster surveys (span- 
ning the redshift range z < 0.5), and simultaneous mea- 
surements of the observable-mass scaling relations, to test 
the consistency of the observed growth rate of clusters with 
that predicted for GR. We allow departures from GR as de- 
scribed by the parameterization Q m _(z) 1 , for which 7 ~ 0.55 
corresponds to GR. We find that the current data are con- 
sistent with GR and give tight constraints on departures 
from it. Assuming constant scatter and purely self-similar 
evolution of the cluster observable-mass scaling relations, 
and a flat ACDM background evolution model, the com- 
bination of XLF, / gas , SNIa, BAO, and CMB data gives 
7(<78 /0.8) 6,8 = 0.55^q '\g. Allowing w to vary, we simultane- 
ously demonstrate consistency with the expansion history of 
ACDM (i.e. w = -1). 

Our analysis employs the method described in Paper I, 
which self-consistently models X-ray cluster survey data and 
follow-up observations to simultaneously constrain both cos- 
mological and scaling relation parameters. This allows us to 
constrain 7 even when allowing for departures from self- 
similar evolution and/or constant scatter in the luminosity- 
mass relation. We caution that the adoption of such a self- 

14 Since, for this figure, we are primarily interested in the ability 
of the CMB data to improve our knowledge of 7 by constraining 
other cosmological parameters, we do not include the additional 
constraining power on 7 from the ISW effect (i.e. we assume that 
for the large scales relevant for the ISW, GR is recovered). This 
has only a small impact on the results (see R09 and Table 1). 



consistent approach to model the survey+follow-up data is 
important, else spuriously tight, and potentially biased, con- 
straints may be obtained. 

Our measurements of the growth of the most X-ray lu- 
minous, most massive galaxy clusters within z < 0.5 con- 
strain the growth of cosmic structure during the epoch of 
cosmic acceleration. It is primarily this redshift dependence 
that provides our constraints on 7, although we have also 
utilized the small contribution from the ISW effect, which 
probes larger scales at z < 2. The inclusion of CMB and 
other data tightens significantly the constraints on the back- 
ground expansion model and other cosmological parameters 
leading to a clear correlation between 7 and as- This high- 
lights the potential for further improvements in 7 with the 
incorporation of precise, independent measurements of as. 
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